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Abstract. We study quantum entanglement in a single-level quantum dot in the linear-response regime. The results 
show, that the maximal quantum value of the conductance 2e^ /h not always match the maximal entanglement. The 
pairwise entanglement between the quantum dot and the nearest atom of the lead is also analyzed by utilizing the 
Wootters formula for charge and spin degrees of freedom separately. The coexistence of zero concurrence and the 
maximal conductance is observed for low values of the dot-lead hybridization. Moreover, the pairwise concurrence 
vanish simultaneously for charge and spin degrees of freedom, when the Kondo resonance is present in the system. The 
values of a Kondo temperature, corresponding to the zero-concurrence boundary, are also provided. 

PACS. 73.63.-b Electronic transport in nanoscale materials and structures - 03.65.Ud Entanglement and quantum 
nonlocality - 03.67.Mn Entanglement production, characterization, and manipulation 



1 Introduction 

' Quantum entanglement, as one of the most intriguing features 

' of quantum mechanics, was extensively studied during the last 
decade, mainly because its nonlocal connotation 1 1 1 is regarded 

I as a valuable resource in quantum communication and informa- 
tion processing 1 2 1 . The question about the relation between the 
entanglement and quantum phase transitions 1 3 1 have been ad- 
dressed recently, for either quantum spin |4 5 6 7| andfermionic 
L8..,9..10J systems, in hope to shed new lights on fundamen- 

I tal problems of condensed matter physics. For example, it was 
shown for spin model |4|, that the entanglement of two neigh- 
boring sites displays a sharp peak either near or at critical point 
where quantum phase transition undergoes. Recently, a class of 
systems with divergent entanglement length away from quan- 

' tum critical point, since the correlation length remains finite, 
was identified |5|. The spin-orbital entanglement analysis was 
also shown to provide a valuable insight to the nature of Mott 

] insulators (6|. In the field of fermionic systems, the local entan- 

■ glement was successfully used to identify quantum phase tran- 
sitions in the extended Hubbard model 1 10|. A separate issue 
concerns using the entanglement as a criterion of quantum co- 
herence LI li when analyzing nonequilibrium dynamics of the 
system with spontaneous symmetry breaking |7 1. 

Here we follow the above ideas, but focus on the physical 
system which undergoes the crossover behavior instead of a 
phase transition: a quantum dot in the Kondo regime. Namely, 
we address the question whether there exist a relation between 
entanglement and conductance for this system? Some earlier 
study mentioned the total entanglement of electronic degrees 
of freedom in the SU{4:) system below the Kondo temperature, 
without determining a qualitative measure of such an entan- 
glement fl2]. In this paper, we consider the SU{2) case, and 



analyze two different definitions of the entanglement between 
quantum dot and the leads: ^raf based on the von Neumann en- 
tropy, and second utilizing the Wootters formula L13J for the 
formation concurrence of two-qubit system. 



2 The model and its numerical solutions 

We study a model of a quantum dot with a single relevant elec- 
tronic level coupled to the left {L) and right (R) metallic elec- 
trodes. The Hamiltonian of the system is 



H = Hl + Vl + Hc + Vr + Hr, 



(1) 



where He models the central region, i?L(fl) describes the left 
(right) lead itself, and Vi(_R) is the coupling between the lead 
and the central region. Namely, we have 



He = edUd + Und^ndi, 



H 



L{R) 



E (4. 



h.c. 



(2) 



V, 



L{R) 



= -Vy c] 



h.c 



Here, Ud — d^d^ is the quantum-dot charge, is the po- 
sition of the molecular level and U is the Coulomb repulsion 
between two electrons. Both Hj^/r) and V[^(r) terms have a 



tight-binding form, with the hopping t and the dot-lead hy- 
bridization parameter V, (cja) creates (destroys) an elec- 
tron with spin cr on site j, the indexes jL{R) denotes termi- 
nal sites of the left (right) electrode. The system is depicted 
schematically in Fig.^ 
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Fig. 1. The Anderson impurity model realized as a double quantum 
dot attached to the leads. The dot is described with the energy level ed 
and the Coulomb interaction U. 

There are many theoretical methods in the existing litera- 
ture, developed to study the electron transport in the presence 
of interaction. In particular, the zero-temperature conductance 
of the quantum dot acting as an Anderson impurity were ob- 
tained within the Bethe ansatz approach 1 14 1 . For the more gen- 
eral situation, one can refer to the Numerical Renormalization 
Group f 15] or to the nonequilibrium Keldysh formalism 1 16|. 
For example, the former approach was succesfuly generalized 
to study a molecule with the electron-phonon coupling 1 17|, 
whereas the latter was adapted for an analysis of the compe- 
tition between the Fano and the Kondo resonance in various 
nanodevices 1 18 1. 

However, since we are interested either in equilibrium trans- 
port properties or in the ground-state quantum entanglement, 
the most useful choice is the variational method recently pro- 
posed by Rejec and Ramsak 1 19 20 1, in which the real-space 
correlation functions are obtained directly. For the system de- 
scribed by the Hamiltonian Q the method converges to the ex- 
act solution II14I . it can also be generalized for multiple quan- 
tum dots 1 19], for the case with a nonzero magnetic field f201, 
or combined with an ab initio wave-function readjustment 121 1 
in the framework of EDABI method 1221 . 




' ' • • • • • ' I 
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Fig. 2. The local entanglement (a) and normalized conductance (6) 
for the system in Fig. Q as a function of the dot energy level td and 
dot-lead hybridization V (changed in steps of O.OSi). 



3 The quantum entanglement 

For the spin s = 1/2 fermionic system, there are four possible 
local states at each site, \v)j = |0)j, | t)j, | \ Ti)i- The di- 
mension of the A^-site system is then 4^ and \vi,V2t ■ ■ ■ tVn) ~ 
rijLi ^re its natural basis vectors. Alternatively, one can 
label the basis vectors by specifying occupation numbers for 
each site and spin \h'i . . . vn) = |"-it ■ • ■ "-att) • ■ ■ "^Ni) , 
with rija = 0, 1. The reduced density matrix for the ground 
state is 

P^a,Ja' =^na,3a'W){^l (3) 

where Tiig-.ja' stands for tracing over all sites and spins except 
the ia and jcr'-th sites. 

3.1 Local entanglement and conductance 

We focus now on the local entanglement |9 1, which exhibits the 
quantum correlations between local state of a selected j-th site 
(e.g. the quantum dot) and the other part of the system (here: 
the leads). For i = j — d and a' = a, the reduced density 
matrix, defined by Eq. (|3}, takes the form 

Pd = u+ 10) (0| + w, IT) (Tl + W2 \i) (il + ITi) (Til , (4) 

where 

= ((1 - "dT)(l - ^dl)) , wi = (ndT(l - ndi)) , 



W2 = ((1 - nd^)ndi) , M- = {ud^Udi) , (5) 

and the averaging is performed for the system ground state. 
Consequently, the corresponding von Neumann entropy Ey (he- 
reinafter called the local entanglement) measures the entangle- 
ment of the states of quantum dot (j — d) with that of the 
remaining — 1 sites, and is given by 

Ey ~ — W+ log2 W+ — W\ log2 W\ — W2 log2 'W2 ~ U_ logj U- . 

(6) 

In Fig. 12 we compare the the local entanglement E^ with 
the conductance calculated from the Rejec-Ramsak two-point 
formula 1191 

where Go = 2e^//i is the conductance quantum, Z\ = 1/Np{ep) 
is the average level spacing at Fermi energy, determined by the 
density of states in an infinite lead piep), E{tt) and £'(0) are 
the ground-state energies of the system with periodic and an- 
tiperiodic boundary conditions, respectively. We found that the 
system size of the order of iV 1000 provides an excellent 
convergence for both the conductance G and the local entan- 
glement Ey (the latter aspect has not been analyzed numer- 
ically before). In particular, the data for N — 1000 cannot 
be distinguished from the ones for N — 2000 in the scale of 
Fig. 12] We also checked that the results are insensitive to the 
number of basis functions composing the Rejec-Ramsak vari- 
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Fig. 3. The local entanglement Ev (a) and spin-magnitude parameter 
0s = (4/3) (Sd) (&) as a function of the dot filling {nd) and dot- 
lead hybridization V . The magnitude of hybridization V goes from 
Q.lt to 0.5t in steps of 0.054. The limiting curves for U/V — and 
U /V = oo are depicted with dashed lines. 



ational wave-function |19|, providing it is ^ 3. When calcu- 
lating correlation functions (|5}, determining the density matrix 
0, one has to choose boundary conditions which minimize the 
ground-state energy for a given system size N: namely, peri- 
odic for N = Ak + 2, and antiperiodic for N — Ak J23J for the 
half-filling l24l . 

Surprisingly, the maximal entanglement between a quan- 
tum dot and leads not always match the maximal conductance 
G = Go = 2e^/h. For small values of the dot-lead hybridiza- 
tion {V < 0.25t for U — t), E„ has a minimum at the parti- 
cle-hole symmetric point — — C//2. It is well known |25| 
that in the limit V'^ /t <C U the Anderson Hamiltonian Q re- 
duces to the symmetric Kondo model with an exchange cou- 
pling p{ep)JK = 8V'^ / ttU \/ At^ — €f- This observation sug- 
gest an important role of the localized moment presence in the 
dot, which strongly affects the entanglement between the quan- 
tum dot and the leads, without observable change to the con- 
ductance. The latter refers to the situation below the Kondo 
temperature Tk oc cxp{—l/ p{ep)JK), above which the con- 
ductance at the particle-hole symmetric point = —U/2 
is depressed and, subsequently, each zero-T Kondo peak dis- 
played on Fig.|2j3 splits into Coulomb-blockade peaks present 
in a finite-T situation |15 16 17 18|. In particular, for U = t 
and V — 0.25<, the exact formula 1.26.1 gives the value of 
Tx/t « 3mK/eV, which seems to be in the experimentally 
accessible range. A finite-T analysis is, however, beyond the 



scope of this paper, since we focus here on T = situation. A 
further discussion of the relation between spin fluctuations and 
the entanglement at the ground state is provided below. 

For the better overview of the system properties we analyze 
them as functions of the dot filling (n^), as displayed in Fig.|3] 
The universal formula for the conductance (not shown) follows 
from the Luttinger theorem 1 27 1 

G = Gosin2(7r(nd)/2), (8) 

whereas Ey evolves gradually from the limit 

i^r" = -(n.)lo& M_(2_(.,))log, (l-^) (9) 
to 

E^-"- = - (l-(nd»log2 (1-M) - M log, 

(10) 

as presented in Fig.|3t. Therefore, in the strong coupling limit 
Ey has two maxima at {ud) =2/3 and 4/3, instead of a sin- 
gle one for {rid) = 1, present in the noninteracting case. The 
major qualitative difference between our results and that ob- 
tained for the extended Hubbard model 1 10| is that Ey behaves 
analytically for any U < oo. This is because the Kondo sys- 
tem, considered here, shows the crossover behavior instead of 
a quantum phase transition present in the Hubbard chain. 

One can also observe, that the entanglement behavior near 
the particle-hole symmetric point (ud) = 1 is determined by 
the magnitude of spin fluctuations, presented in Fig.|3j5- As a 
measure of such fluctuations, we choose a parameter I 28J 

Os^^iSl) ^{nd)-2{nd^ndi), (11) 
which obeys the inequality 

(n,) (12) 

where the lower and the upper limit refers to the U — and 
U = oo case, respectively. In particular, for (n^) = 1 the spin- 
fluctuation parameter varies from 6*5 = 1/2 for free fermions 
to 6*5 = 1 for the localized spin-1/2. The charge fluctua- 
tions are determined by 0s as Varjnrf} = (n^) — (ud)'^ — 
{nd){2 — {ud}) — Os, so for (n^) = 1 and [/ = 00 we ob- 
tain Var{nd} = 0. The vanishing of charge fluctuations and 
the value of the spin-square (S^) = 3/4 allows one to con- 
sider the localized spin 1/2, a presence of which governs the 
ground-state properties at the strong-coupling limit. The den- 
sity matrix takes thus the form pd = (IT) (T| + |i> (i|)/2, 
which brought us to the value of the local entanglement Ey — 1 
at the particle-hole symmetric point. In contrast, for the non- 
interacting system all coefficients of the density matrix are 
equal to 1/4 and the local entanglement reaches its maximal 
value Ey = 2 (for (ud) = 1). The entanglement drop with the 
increasing coupling near the particle-hole symmetric point can 
therefore be explained as an effect of the formation of a local- 
ized moment inside the dot. 

The coiTespondence between inequalities il2\ and the lim- 
its defined by Eqs. (|9} and (I10> become straightforward when 
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expressing the coefficients of the density matrix as func- 
tions of {nd) and 0s, what leads to the local entanglement 



{rid) 



Os 



9s 



log; 



f ind) 



log2 



Os 



2 - {ud) - Os 



Os log; 



Os 



(13) 



Eq. (I13> with the limits given by ( I12t relates the entanglement 
between the quantum dot and the leads to the local moment 
formation inside the dot. It also express the local entanglement 
Ey in terms of measurable quantities: the dot occupation (ud) 
and the spin-square magnitude (S^) contained in the parame- 
ter 6*5 il 1> . In contrast, the spin fluctuations are absent in Eq. 
(|8j for the conductance G, which is fully determined by the 
dot filling (ud) ■ One can note Eq. il3\ is model-independent, 
providing we consider the lattice system with one orbital per 
site. Thus, for the system with quantum phase transition, such 
as that considered by Gu et al. 1 10|, the nonanalytical behavior 
of the local entanglement Ey is equivalent to the nonanalytical 
behavior of the spin-magnitude parameter Os- 



3.2 The fermionic concurrence 

We consider here the entanglement of two qubits, one associ- 
ated with the electron localized on a quantum dot and other 
with the nearest one placed in a lead. The physical realization 
of individual qubits may, in principle, employ charge or spin 
degrees of freedom of the system in Fig.Q 

The reduced density matrix Q for the pair of electrons with 
equal spins (say a = a' =1), one localized on a quantum dot 
and other on a nearest lead atom (i — d, j = iL(R))^ can be 
written as 1291 



■'I I J I 





wl z"" Q 
{z*Y 
ui 



(14) 



where 



((1 



)(1 



(1 



W2 = ((1 - "-iT)"-jT)> = (c]t'^'T)' = (15) 

The upper index c stands to denote, that the density matrix MA\ 
refers to the charge degrees of freedom, since the spin direction 
is arbitrarily chosen for both particles. 

We use now the concurrence C as a measure of the entan- 
glement for such a two-qubit system. The closed-form expres- 
sion, derived by Wootters II13I . reads 



C 



0,\/^i - 



A., - J\ 



(16) 



The Ai's are the eigenvalues of the matrix product p ■ (cr. 



V 



(Jj)p*{(Tf (8) crj), and Ai ^ A2 ^ A3 ^ A4. Since there exists a 
monotonous relation between the concurrence C and the entan- 
glement of formation Ef — ~x log2 x — {1 — x) log2(l — x). 
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Fig. 4. The charge (a) and spin (b) pairwise concurrence for the sys- 
tem containing one qubit localized on quantum dot and other on the 
nearest atom of the lead. 



where x = 1/2 + Vl — C^/2 1 13 1, C is widely used instead of 
Ef m the literature. 

For the density matrix pi]j] J14t the corresponding con- 
currence can be calculated from Eq. ( I16> as 



Ci^^j^ — 2 max 



|0, {c\^c.j^ 

Kt"jT> (1 - ("»T> - Kt> + ("»T"jT>) 

(17) 

Hereinafter, we call Ci|.j| a charge concurrence, since it is re- 
lated to the charge degrees of freedom. 

Alternatively, one can consider the full two-site density 
matrix pij = Tr^ jtZ/) (tf^l (with i = d, and j = ji(_R,) again) 
and project out all the states except from these coiTesponding 
to fii = Uj = 1. The resultant 4x4 density matrix Ps^.s, ~ 
Psi.sj /TiPsi.sj describes the entanglement accessible by spin 
manipulation with a particle-conservation constrain |30|. The 
matrix ps- has a general structure of Pi^j^ given by Eq. il4\ . 
with the elements (I15> replaced by 

u+ = (?^^T(1 - "a)»^iT(l - "ji)) ' 



/ = ((1 - njT)"a(l - '^jT)"ji> ' 
«i = {nz^i^ - "a)(l - '^jT)"-ji) ' 
'^2 = ((1 - n,T;)niinji {l - n^i)) , 



(18) 
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Table 1. The critical values of hybridization Vc on the zero- 
concurrence boundary for = —U/2, the corresponding exchange 
coupling p{0)Jk = 4V'^ /ivUt, and the Kondo temperature L26J . 



U/t 


14 


p(0)Jk 


TK/t 

[mK/eV] 


0.5 


0.23 


0.135 


360 


1.0 


0.26 


0.086 


8.6 


2.0 


0.31 


0.061 


0.13 



that the coexistence of zero concurrence and maximal conduc- 
tance can be observed for U ^ tif the relative temperature T/t 
is of the order of a few mK/eV. Another interesting feature of 
these results is the universal (interaction independent) behavior 
of the maximal dot filling (ud) = ("-d)max, for which the en- 
tanglement C ^ 0, with V/t (cf. Fig.|5j'). This observation 
can be rationalized by using Eq. (I17> and putting Ciui — 0. 
Then, in the perturbative limit F <C t, we obtain 

(nd)max«23/2|(4c,)|«25/2(FA)2, 

up to the quadratic terms. The agreement with the numerical 
data is perfect for V/t ^ 0.1. 



Fig. 5. The values of dot energy level (a) and average occupa- 
tion (ud) (b) corresponding to zero concurrence, as a function of hy- 
bridization V and Coulomb interaction U (specified for each dataset). 
The perturbative limit V <^tis also shown. 



The label s indicates that we are working know with the spin 
degrees of freedom, as charges of the i and j sites are chosen. 
The concurrence, obtained by applying the definition il6\ to 
the density matrix Psi.s -, is called a spin concurrence. 

The charge and spin concurrence is shown in Fig.l^as a 
function of the dot energy level e^. Again, we observe an excel- 
lent convergence of both the studied quantities for the system 
size of the order of iV 1000. Although the charge and spin 
concurrence are, in principle, two different physical quantities, 
they reach the limit C — Q simultaneously for all the analyzed 
values of V and U . Therefore, we can conclude that below the 
critical value of the hybridization V < Vc{U) and in the Kondo 
regime, the qubit localized on the quantum dot is not entangled 
with other placed on the top of the lead for neither charge nor 
spin degrees of freedom. This is because of large quantum fluc- 
tuations in each lead (modeled as a noninteracting Fermi gas), 
which therefore destroy the entanglement for a weak dot-lead 
coupling. 

The values of and {ud), corresponding to C = are de- 
picted in Fig.|5l To complement the analysis we also provide, in 
Table 1, the values of the symmetric Kondo temperature |26|, 
corresponding to the critical hybridization Vc{U), at which the 
concurrence vanish in the ground state. Although the numeri- 
cal results presented in Fig.|5lrefers to T = 0, it is clear that in 
the finite-T situation the conductance G Gq^A the particle- 
hole symmetric point if T <C Tk, and that the concurrence C 
decrease with T due to dephasing. Therefore we can conclude. 



4 Summary 

We analyzed the local entanglement between the quantum dot 
and the leads as a function of the dot energy level e^, the dot- 
lead hybridization V and the intra-dot Coulomb repulsion U . 
The measure of this entanglement, the von Neumann entropy 
Ev, evolves gradually from the weak-coupling limit, in which 
the maximal _E„ match the maximal quantum value of the con- 
ductance G = Go = 2e^ /h, to the strong-coupling situation, 
where maximal G corresponds to the local minimum of Ey. 
This behavior was explained in terms of local moment forma- 
tion inside the dot, which took place when the charge transport 
is dominated by the Kondo effect. 

Finally, we defined the pairwise concurrence, measuring 
the entanglement between a pair of qubits: one localized on 
the dot and other on the nearest atom of the lead, for charge 
and spin degrees of freedom separately. Both quantities van- 
ish simultaneously in the Kondo-resonance range, where the 
weakly-entangled system show the maximal conductance. We 
predict the latter to be observable at the temperature range of 
T/W ^ 1 mK/eV (where W = 4t is the lead bandwidth), 
which seems to be accessible within the present nanoscale ex- 
perimental techniques. The universal dependence of the maxi- 
mal dot filling, above which the concurrence vanish, {nd)^^^^^ ~ 
25/2(T//i)2 for V/t < 1, was also identified. 
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